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Abstract 

A supersymmetric standard model with heavier scalar particles is very interest- 
ing from various viewpoints, especially Higgs properties. If the scalar mass scale is 
O(10 - 10 4 ) TeV, the standard model-like Higgs with mass around 125 GeV, which 
is implied by the recent LHC experiments, is predicted. However this scenario is 
difficult to be directly tested with collider experiments. In this paper, we propose a 
test of this scenario by using observations of primordial gravitational waves (GWs). 
The future GW experiments such as DECIGO and BBO can probe the scalar mass 
around 0(1O 3 — 10 4 ) TeV, which is preferred from the Higgs mass about 125 GeV, 
if the primordial GWs have large amplitude. 



1 Introduction 



A supersymmetric (SUSY) standard model (SSM) is a very appealing model beyond the 
standard model (SM). An important prediction of the SSM is lightness of SM-like Higgs 
particle. The recent observation by ATLAS and CMS collaborations suggests the existence 
of the Higgs particle with mass around 125 GeV [HE]. To realize the Higgs mass about 
125 GeV, the SSM can be strongly constrained: heavy scalar top quark and/or very large 
A-term of the scalar tops are required in the case of the minimal SSM (MSSM) [31 HJ E]. 

Among them, heavy scalar scenarios are interesting from various viewpoints. With 
the small A-term of the stops, stop mass around O(10 — 10 4 ) TeV can successfully explain 
the Higgs mass about 125 GeV [5j El [7]. In addition to the Higgs mass, this scenario is 
successful in the SUSY flavor and CP problem and/or the cosmological gravitino problem 
[8] as discussed in the context of the split SUSY pi [101 E] or the anomaly mediated 
SUSY breaking (AMSB) models [121 H3]- Motivated with such a scenario with heavy 
scalar mass, phenomenological consequences are discussed in Refs. [HE]- 

However, this scenario is very difficult to be directly tested with collider experiments. 
This is because the scalar particles are too heavy to be produced at colliders, even if 
the gauginos are light enough to be produced at colliders. However once the gluino-like 
particle is discovered, there is a chance to probe the scalar sector. The gluino decay is very 
sensitive to the scalar sectors. For example, if the scalar mass scale is larger than about 
10 3-4 TeV, the gluino has a sizable lifetime. Therefore, an R-hadron or a displaced vertex 
of a long-lived gluino can be a probe of such a heavy scalar mass [HI US] - Although this is 
an interesting possibility, it is still indirect evidence for the heavy scalar SUSY scenario. 
Thus, it is very important to investigate the independent way to probe the scalar sector. 

In this paper, we discuss testability of this scenario with the next generation projects 
of gravitational wave (GW) experiments as Deci-hertz Interferometer Gravitational Wave 
Observatory (DECIGO) pU OH [18] and Big Bang Observer (BBO) pj5]. A main target 
of the next generation GW projects is the primordial GWs generated in the inflationary 
era [20]. The observations of the primordial GWs can provide information not only on 
inflation but also on the thermal history of the early Universe [211 [22]. In the case of 
the heavy scalar top mass scenario, it is natural to expect that the other SSM scalar 
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particles have masses of the same order of magnitude. Therefore, the sudden change of 
relativistic degrees of freedom as much as ~ 100 is expected at temperature around the 
scalar mass, which leads to a step at the corresponding frequency in the amplitude of the 
GWs [221 123]. The relevant energy range O(10 2 — 10 4 ) TeV corresponds to the frequency 
band of DECIGO/BBO, 0(0.001-0.1) Hz. We then investigate whether DECIGO/BBO 
can detect a step in their band caused by the change in the relativistic degrees of freedom 
using the correlation analysis techniques [24]. By measuring frequency dependence of 
the primordial GWs in detail, we can obtain information of very high-mass-scale particle 
beyond the reach of the collider experiments, which gives a good test of the SUSY models. 

2 Heavy scalar scenario and cosmology 

We describe here the Higgs properties and cosmological consequences of our setup with 
high SUSY breaking scale. In the low-energy below the heavier SSM scalar masses, there 
is a SM-like Higgs scalar doublet H. The relevant Higgs potential is similar to the SM 
one and written as 

V(H) = ±(Htf-v*y, (1) 

where v — 174 GeV is the Higgs VEV and the Higgs mass is provided as rrih = \/2\v. At 
tree-level in the MSSM, the A comes from D-term potential and is related to the gauge 
couplings as, 

A = ^(^ 2 2 + ^)cos 2 (2/3). (2) 

Hence, the value of the Higgs mass cannot be arbitrarily large. Instead, the Higgs mass 
is bounded above as, 

m h = m z cos(2/3) < 91 GeV, (3) 

at tree-level where f3 is defined by tan/3 = v u /v<i- This value is quite small compared to 
plausible Higgs mass ~ 125 GeV. A large value of the Higgs mass is realized by including 
contributions from the radiative corrections, Am^, which are controlled with the SUSY 
breaking scale [3]. Although the detailed value depends on SUSY breaking parameters 
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and tan /3, in the case of small A-term of the stops, which is plausible in the context of 
the split SUSY and AMSB scenarios, O(10 — 10 4 ) TeV stop mass provides appropriate 
Higgs mass. In such a case, it is natural to expect that the masses of MSSM scalar other 
than stop are also heavy, m sca i ar ~ m stop . For example, in the case of the AMSB models, 
the scalar masses are controlled with the gravitino mass m 3 / 2 and, hence, all the scalar 
masses are likely as large as the stop mass. Therefore, there are expected to be nearly a 
hundred of particles other than the SM above the scalar mass scale, m sca i ar . Hereafter, 
we assume that all the MSSM scalar particles have universal scalar mass mo and that 
the /x-term is also large: fi = m . These assumptions are just for simplicity. Even if we 
consider the variation in the scalar masses, the results presented below do not change so 
much. 

Then let us consider the cosmological history in this scenario. In this paper, we focus 
on a fact that the effective number of relativistic degrees of freedom changes at the scalar 
mass scale, T ~ m . We define the effective number of relativistic degrees of freedom 
and g* s as, 



P =^T\ (4) 
. = ^tflr, (5) 

where p and s are the energy density and the entropy density, respectively. The evolution 
of g* and g* s are determined by the MSSM mass spectrum. In Fig. [TJ we have shown an 
example of g* and g* s when m gaug i no = 1 TeV and m = 10 4 TeV. From the figure, we can 
see that the relativistic degrees of freedom increase by 100 at temperature corresponding 
to the scalar mass scale, T ~ 10 4 TeV. 



The g* change discussed here causes features in the GW spectrum [22} [23] . In the 
next section, we investigate the detectability of this g* change by the next generation GW 
experiments as DECIGO pH [18] and BBO p2]. 
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Figure 1: Evolution of and g* s as functions of temperature T. We have set m 
TeV and m = 10 4 TeV. 
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3 Detecting change in the gravitational wave back- 
ground spectrum 

3.1 Energy spectrum of gravitational waves 

We first describe the feature in the GW spectrum expected in this scenario. The amplitude 
of the stochastic GWs is conventionally provided in terms of the relative spectral energy 
density defined as, 

ugwC/) - -p^- (6) 

Pcdlog/ 

Assuming the primordial tensor power spectrum is almost flat in the relevant frequency 
band, the amplitudes for modes crossing the horizon before and after the change are 
related as [22], 

^GW (/after) = ( ^*' aftCr j f2 G W (/before) , (7) 
\9*, before / 

assuming the horizon crossing occurs in the radiation-dominated era since the GW en- 
ergy density pew evolves as oc a~ A after the horizon crossing and the critical density p c as 
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oc g* a -4 , respectively. Here, we have assumed g %s = g* because their difference is irrel- 
evant to the present analysis. Then, we expect about 20% change in the GW amplitude 
in this scenario. 

The step appears at frequencies corresponding to modes that cross the horizon during 
<7* change. For modes crossing the horizon in the radiation-dominated era, the frequency 
is related to the horizon- crossing temperature T c as |22j . 

/ = 2.7 x 10" 5 Hz ( -^—) (^-)\ (8) 
J \1 TeVy U28/ 1 ; 

where 128 corresponds to the approximate value of g* of the SM particles and the gauginos 
(see Fig. [I]). Since the relevant energy range O(10 2 — 10 4 ) TeV corresponds to the 
frequency band of DECIGO/BBO, 0(0.001 - 0.1) Hz, DECIGO and BBO have a potential 
to detect the g* change expected in this scenario. 

In the slow-roll inflation paradigm, the GW spectrum ([6]) is approximated to be flat 
within the DECIGO/BBO band. Then, the GW spectrum can be expressed as, 

i 

3 

^GW,0) (9) 

where g*(T c ; mo) is the number of relativistic degrees of freedom provided in the previous 
section for the scalar mass m Q and T c is related to the frequency as Eq.®. Though the 
spectrum could have a tilt within the DECIGO/BBO band, we can show that our results 
are not affected much even if we consider the tilt of the spectrum given by the standard 
inflation consistency relation between the GW amplitude and the spectral index |25j . 

The GW spectrum OH]) is characterized by three parameters, which we have denoted as 
6: the amplitude ficw.O) scalar mass m , and the value of g* at high temperature <?*,before, 
which respectively correspond to amplitude in the low-frequency region, frequency at 
the step, and amplitude in the high-frequency region. For convenience, we employ a 
normalization related to the tensor-to-scalar ratio as, 

fi GW ,o = 4.4 x 10- 15 r GW , (10) 

though the analysis is independent of the CMB observations. The relation between r G w 
and v CMB; the tensor-to-scalar ratio on the CMB scale, depends on the model of inflation 



g*(T c ;m ) 
128 
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(the primordial spectrum for tensor modes) [26]. The value of row is typically (9(1 — 10)% 
less than tcmb for relevant cases though some models predict larger values [27]. The 
relation also depends on the thermal history of the Universe. In this scenario, for example, 
additional entropy production is expected from decays of some meta-stable particles such 
as a gravitino. However, this contribution hardly affect the relation to an accuracy of the 
analysis in the typical SUSY models. Since r cmb is constrained as r CMB < 0.36 at the 
95% level of confidence [28], row is expected to be 0(0.1) at most. We also introduce the 

Change in as Ag* = #*,before - 0*,after = 0*,before - 128. 

In the remaining part, we investigate detectability of the step induced by the g* change 
with DECIGO and BBO for scalar mass mo in the range (9(10 — 10 4 ) TeV. 

3.2 Detectability of the step 

To detect a weak signal as the primordial GWs, we should perform a correlation analysis 
between two detectors. Here, we investigate whether DECIGO/BBO can detect the step 
by using the correlation analysis techniques based on Ref. |24j . 

The output for each detector in the Fourier space £>i(f) (i = 1,2) are written in the 
form, 

&(/) = W) + *(/), (ii) 

where hi(f) is the noise and Sj(/) is the contribution from GWs. Dividing the frequency 
band into segments = [fa — 4r, + 4r] (i = l,...,iV b ), we consider the following 
combination of the outputs for two detectors, 

*™ s ^Kf ® L if mmf) ' (12) 

where T Q b s is the observation period and Hq = 73.5 km/s/Mpc is the Hubble parameter 
at the present time [28] . Taking Af sufficiently small, we can evaluate the GW spectrum 
from the expectation value of Ogw,;, 

(n GW ,i) ^ n G w(/i), (13) 
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with the variance, 
1 



Afi 



GW. 



.2 _ 



T obs A/ l 
Afi GW (/, A2 



2fi G w(/,) 2 + 2cl ( 8 ^£j f?n G M)s n ( fi ) + 4 2 /« s^) 2 



(14) 



where S n (f) is the square spectral noise density. Here we have introduced a factor c n 
to see how improvements in the strain sensitivity, a/S^, change the detectability. If we 
take the weak-signal limit, which is appropriate for near- future detectors, the expression 
is simplified as, 

^^S'w (15) 

In Fig. [21 we have shown the GW spectrum with m = 10 4 TeV and Ag* = 100 and 
the noise curves for DECIGO, BBO and ultimate-DECIGO p2J [31]. Here, we depicted 
(A///) 1 / 2 Afi GW defined in Eq. (JUD as the noise and we set T ^ s = 1 yr. We have used 
fitting formulae presented in Ref. [29] for estimating the square spectral noise density. In 
the case of ultimate-DECIGO, the noise (JHJ) is dominated by the contributions from GW 
itself at lower frequencies. 

Detectability of the step is determined by whether the spectral shapes with different 
parameters, = {row, m o, Ag*} and Of = {rcw,/, m o,f, Ag*j}, can be distinguished from 
the data ([12]) . To quantify the detectability, we introduce the following quantity, 

m s '> = Af JJ f ^kMff ' (16) 

where is the number of correlation signal, which we set = 2.0 Intuitively speaking, 
T represents how Ogw(/;^/) an d Ogw(/; ^) differ from each other relative to the noise 
Af2 G w- 111 short, the quantity J 7 characterizes the expectation that the true parameters 
Of are confused with in the future experiments. We use the value of J 7 to estimate the 



confidence level (C.L.) in the future GW experiments. See Appendix A for a justification. 



In Eq. (j!6p . we have introduced a lower cut-off frequency / cut to avoid the confusion 



1 From Eq. (fT()]) . we can see that (A/ / f) 1 ^ 2 AQq-w provides effective noise in the present anaiysis for 
Slow per log /, which we have plotted in Fig. [2] 
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Figure 2: The GW spectrum for m = 10 4 TeV (black line) and noise curves for DECIGO 
(blue), BBO (red) and ultimate-DECIGO (green). Here we set T obs = 1 yr. Note that 
the noise (|14|) is determined by the contributions from GWs at lower frequencies for 
ultimate-DECIGO. 

noise from the cosmological white dwarf (WD) binaries [30]. Strictly speaking, we should 
also consider contributions from the neutron star binaries. Here, we simply neglect them 
because they can be subtracted from the data for the sensitivities required for detecting 
the g* change as we see below [29J. The noise introduces a cut-off in the high frequency 
region, / max , to Eq. (TTS]) . though we set / max = 100 Hz in the estimation. Note that T 
does not depend on the bin size Af. Thus, the value of Af is irrelevant to the analysis if 
it is sufficiently small. We discuss the detectability by evaluating J 7 for Ag*j = 100 and 



In Fig. [3j we have depicted J 7 minimized with respect to rcw^cl J^mm, for Ag*j = 100 
and Ag* = 0. Here, J 7 is minimized to take into account that degeneracy with rcw 
increases estimation error. The curves are presented for various values of the improvement 
factor c n = 1, 1/2, 1/3. In the estimation, we have set the other parameters as ?"gw,/ = 0.1, 
T bs = 5 yr, and f cut = 0.1 Hz. Note that the values of .F m in for another parameter choice 
can be easily obtained by scaling J-" m i n as oc T b s ?"Q W jC~ 4 in the weak-signal limit. We 

2 Note that T m i n is independent of m because the GW spectrum is flat for Ag* = 0. Therefore, J^nin 
is given by a function of m/ for Ag*j = 100 and Ag, = 0. 



Ag* = 0. 
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have also shown the curves for different lower cut-off frequencies / cut in Fig. H] to see 
the effect of the WD noise for the estimation. We assume row,/ = 0.1 and T D b s = 5 yr 
and set c n = 1/3 for DECIGO and BBO. To see a detectability in an ideal case, we have 
also depicted J-" m i n for ultimate-DECIGO. In this case, the variance (1141) is dominated 
by the contributions from GWs and the detector noise S n only determines an effective 
cut-off in the high frequency region. Hence, the values of F m i n are approximately constant 
for different values of row./ as far as the contributions from the detector noise can be 
neglected. 

From the figures, we can see that BBO can test our scenario with m = O(10 3 — 
10 4 ) TeV if row is large enough, though testability in the low mass region depends on the 
lower cut-off frequency associated with the WD noise. Though sensitivity of the standard 
DECIGO is not enough to detect the g„ change [23], DECIGO can be comparable with 
BBO if we employ an improved design with three times better sensitivity. A degree of 
improvement in sensitivity is important to extend the testable mass range. A large region 
in the mass range O(10 3 — 10 4 ) TeV can be covered by the interferometers as BBO-grand 
(c n ~ 1/3) or ultimate-DECIGO, though the WD noise should be removed to test the 
lower mass region mo < O(10 3 ) TeV. 




io 3 io 4 io 5 io 3 io 4 io 5 

m„j [TeV] mo,/ [TeV] 



(a) DECIGO (b) BBO 

Figure 3: The measure T minimized with respect to tqw f° r DECIGO (left) and BBO 
(right) for the improvement factor c n = 1, 1/2, 1/3 from bottom to top. Here, we have 
set row,/ = 0.1, T D b s = 5 yr, and / cut = 0.1 Hz for both detectors. The horizontal line 
indicates = 2. 
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(a) DECIGO with c n = 1/3 




m ,/ [TeV] m 0> / [TeV] 

(b) BBO with c n = 1/3 (c) Ultimate-DECIGO 

Figure 4: The measure J 7 minimized with respect to row for the lower cut-off frequency 
/cut = Hz, 0.1 Hz, 0.2 Hz. Here, we have set row,/ — 0.1, c n = 1/3, and T obs = 5 yr. The 
horizontal line indicates J-mm = 2. We have also depicted J-mm for ultimate-DECIGO. 
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In Fig. [51 we have also depicted expected errors in (mo, Ag*) for BBO with c n = 1/3 
from the value of JF-mm-, assuming (r GWi /, m 0) /, Ag*j) = (0.1, 10 4 TeV/3 x 10 4 TeV, 100). 
The other parameters have been set as T ^ s = 5 yr and / cut = 0.1 Hz. Note again that the 
value of J-min for another parameter choice is obtained by scaling J-" m i n as oc T b s ?"Q W jC~ 4 
in the weak signal limit. From the figures, we can see the interferometers as BBO-grand 
have a potential to determine the scalar mass scale, though precise value of error depends 
on the MSSM spectrum. The parameter degeneracies are observed for mo,/ = 10 TeV 
because the WD noise covers the spectrum in the lower frequency region. 




mo/mo j [TeV] rn a /m j [TeV] 



Figure 5: Expected error in (m ,Ag*) with (row,/) m o,f, At?*,/) — (0.1, 10 4 TeV/3 x 
10 4 TeV, 100) for BBO with c n = 1/3 (left: 10 4 TeV, right: 3 x 10 4 TeV). Here, 7 is 
minimized with respect to r^w to take into account that degeneracy with row increases 
estimation error for (mo, Ag*). Here, we have set T Q ^ S = 5 yr and / cut = 0.1 Hz. 



4 Conclusion and Discussion 

In this paper, we have discussed a GW test of the SSMs with heavy scalar particles 
motivated with the recent LHC Higgs search. We have quantitatively shown that the 
future GW experiments as DECIGO and BBO have a potential to test the heavy scalar 
scenario, m = O(10 3 — 10 4 ) TeV, and determine the scalar mass scale. Even in the case 
that the GW amplitude is large, a degree of improvement in sensitivity is necessary to 
cover a large region in the mass range O(10 3 — 10 4 ) TeV by DECIGO and BBO. We 
have seen that it can be covered by using BBO-grand or ultimate-DECIGO. Our analysis 
3 See | Appendix A| for the precise meaning of the confidence level depicted here. 
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depends on the assumption for the spectral shape of the primordial GW spectrum in the 
DECIGO/BBO band as is naturally expected. For example, if the reheating temperature 
to evade the gravitino problem is too low, the spectrum is expected to be deformed in the 
DECIGO/BBO band [32]. However, in the case of very high reheating temperature sce- 
nario, such as the thermal leptogenesis scenario, we can safely make the analysis assuming 
the flat spectrum. In addition, the spectrum could have a tilt in the DECIGO/BBO band. 
Although it can be shown that our results do not change much under the standard in- 
flation consistency relation, introduction of the tilt could increase error in the estimation 
much if no prior knowledge on the tilt is assumed. In that case, better sensitivity in the 
high frequency region will be needed to break the degeneracy of the models. 
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Appendix A The meaning of T 



Here, we clarify the meaning of J 7 given by Eq. ( |T6|) . Provided data Oqw^ the parameters 
9 can be estimated by maximizing a likelihood function provided as, 



GWj 



J\f 
1 

A7 



exp 



exp 



~2~^ 



AO 2 



2* Af 



df 



(n GW (f)-n GW (f;9)y 



(17) 



where M is a normalization constant. □ Here, we have assumed that each segment 
contains so large number of Fourier modes that ficw,; can be considered to be Gaussian 
distributed. This assumption can be compatible with small bin size because the frequency 
resolution, T^ s ~ 0(1O~ 8 ) Hz, is much smaller than the relevant frequency, 0(0.1) Hz. 



A factor 2 appears in Eq. (fT7f because Fourier modes are denned in (— oo, oo). 
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Then, the estimated value of parameters 9, maximum likelihood estimators (MLE), are 
obtained by maximizing the likelihood function <\17L or minimizing the quantity, 



2 _2iV d f (n GW (f)-n GW (f;9)) 2 



A/ L Afi GW (/) 



i. mm ' 



and the error (confidence level) in the parameters by seeing the value of A% 2 = \ 2 ~ xl 

Now we consider the expectation that parameters 9 are judged to be consistent by 
future experiments when true values of parameters are 9f. (^gw(/)) = ^Gw(/j To 
estimate the accuracy of the estimation, we should consider not only the error expected 
in the experiments but also deviation of MLE from the true values, which depends on 
the data obtained in the future. This uncertainty of MLE can be included by taking the 
ensemble average over realizations of the data Ogw(/) with (Ogw(/)) — ^gw(/; #/) : 



(£(0;ft G w)) oc exp 



±.*k [ d f (°gw(/; 9 f ) - n GW (/; 9)) 



2 A/ J f A0 OW (/) 



(19) 



The exponent is just —JF/2. The quantity J 7 measures the compatibility of the future 
data with the model where values of parameters are 9 when the true values are 9f. 

We can relate J 77 to the confidence level expected in the future GW experiments as 
e~ T l 2 = 2 Nd / 2 (e^min/ 2 ' e ~ Ax ^ 'I 2 ) . Here, the average of A% 2 is taken with a weight e~ x ^ 2 , 
the maximum value of the likelihood function, which measures the goodness-of-fit of 
our model with a realization of the data. Hence, J 7 approximates A% 2 averaged over 
realizations except those less compatible with our model, for which the estimation of 
parameters is less reliable. In Fig. [61 we have depicted the result of Monte Carlo simulation 
for parameter regions where parameters are expected to be in the 99% C.L. contour with 
a given probability, assuming (r G w,/i m o,/ ; ^9*,f) = (0.1, 10 4 TeV/3 x 10 4 TeV, 100) for 
BBO with c n = 1/3, T Q b s = 5 yr, and / cut = 0.1 Hz. From the figure, we can see that the 
"99% C.L. contour" depicted by using the value of J 7 covers almost all of parameters that 
can have the 99% C.L. in the future analysis; the probability that parameters outside the 
contour have 99% C.L. is less than ~ 5%. Therefore, we use the value of J 77 , which can be 
estimated analytically, as a measure for the expected error of the parameter estimation 
in the future GW experiments. 
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Figure 6: The result of Monte Carlo simulation for parameter regions where parameters 
are expected to be in the 99% C.L. contour in the future GW experiment. The color 
panel shows the probability that parameters are contained in 99% C.L. contour. The 
dashed and red line represents "99% C.L. contour" depicted by using the value of J 7 and 
the solid and black lines the contours with a constant probability. To depict the figure, 
we have used BBO with c n = 1/3, T Q b s = 5 yr, f cnt = 0.1 Hz, and (row,/, m o,/, ^*,/) = 
(0.1, 10 4 TeV/3 x 10 4 TeV, 100). (left: 10 4 TeV, right: 3 x 10 4 TeV) 
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